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Abstract 

In this paper we review some results on longitudinal beam dynamics obtained in the 
framework of the Thermal Wave Model (TWM). In this model, which has recently shown the 
capability to describe both longitudinal and transverse dynamics of charged particle beams, 
the beam dynamics is ruled by Schrodinger-like equations for the beam- wave-functions, whose 
squared modulus is proportional to the beam density profile. Remarkably, the role of the 
Planck constant is played by a diffractive constant e, the emittance, which has a thermal 
nature. 


1 Introduction 

Recently, on pure basis of analogy with other similar subjects of physics, a new technique to derive 
the equation of motion for a thermal system, like a charged particle beam at finite temperature, 
which is able to take into account the collective behaviour of the ensemble has been obtained 

[l]-[5]. 

The starting point of this technique, the Thermal Wave Quantization (TWQ), are the equations 
of motion of the considered system, in the so called single-particle approximation, from which is 
possible to obtain the single-particle hamiltonian of the system. At this point, the formal anal- 
ogy showed in the case of the transverse dynamics for relativistic charged particle beams, with 
the electromagnetic optics in paraxial approximation and with the two dimensional nonrelativis- 
tic quantum mechanics [1], suggests to replace the single -particle hamiltonian , with a differential 
operator, and the hamilton-equations with a Schrodinger-like equation, in which coordinate and 
particle momentum are replaced by a beam-wave-function. 

This technique, applied to the longitudinal and transverse beam dynamics, has led to the for- 
mulation of the Thermal Wave Model (TWM) for relativistic charged particle beam propagation, 
which represents a useful quantum-like description of the total beam optics [l]-[5]. This model has 
already been successfully applied for estimating the effects of the aberrations in linear colliders 
[3], [5], as well as for describing nonlinear beam-plasma interaction [2], and nonlinear longitudinal 
dynamics in circular accelerating machines [4], 
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2 TWM for longitudinal dynamics 


Let us consider a single relativistic particle of electric charge q, within a stationary bunch, travelling 
with longitudinal velocity (5 c (/? ~ 1) in a circular accelerating machine of radius Rq = cT 0 / 2n (T 0 
being the revolution period). Its longitudinal motion is described, neglecting radiation damping 
and quantum excitations, by a pair of equations which, defining s = ct (t being the time), can be 
put in the following dimensionless form [6] 


dx 

ds 

ctP 

ds 


A E 

V E 0 ~ 
q AV 
cTqEq 


rfP 




( 1 ) 
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where x is the longitudinal particle coordinate and V = A E/ Eq is the dimensionless longitudinal 
energy variation, after a turn in the ring. Note that x (—'kRq < x < ttRq) and V are both evaluated 
with respect to the synchronous particle (A E = 0), and E 0 is the synchronous particle energy. 
The quantity AV represents the total voltage variation after a turn and it takes into account 
the interactions of the particles with the surrounding medium (RF-cavities, pipe, kickers, etc.). 
Consequently, the equations (1) and (2) describe the longitudinal bunch dynamics on time scale 
t » T 0 . Furthermore, in (1) the parameter 77 is defined as rj = (Aw/w 0 )/(AF/Fo) ( w o = c/Ro 
and Au; being the frequency shift with respect to u»o)- By defining the momentum compaction 
a = (AR/Rq)/(AE/Eq), where A R is the orbit radius variation with respect to Ro , it can be 
easily proved that rj = I/7 2 - a. From (1) and (2) we can easily write the following dimensionless 
single-particle hamiltonian 

H = i T)P 2 + U , (3) 

where 1 

u = / gAV dx' . (4) 

cTqEq Jo 

In order to find an equation which describes the longitudinal evolution of the beam, taking into 
account its thermal spreading (longitudinal emittance) while it interacts with the surrounding 
medium (potential well and wake fields), we follow the quantum- analogy, which suggests to use in 
(3) the following correspondence rules 

d d 

V -* V = -itLTT and H —> H = ic L — , (5) 

ox os 

where e L is the longitudinal beam emittance. Thus, by considering (3) and (5), the following 
Schrodinger-like equation for the beam wave function (bwf) 'k can be assumed 

ie L — = H<H , ( 6 ) 

os 


where H = %V 2 + U. Consequently, (6) becomes 



elrjdH 
2 dx 2 


+ (/'k 


(7) 
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Note that (7) describes the longitudinal beam dynamics in terms of the bwf which we as- 
sume to be related to the longitudinal density number \(x,s) through the following relation: 
A(x,s) - A 0 |vf'(x,s)| , where Ao = N/R 0 (N being the total number of particle in a bunch). 
According to the previous definitions, |^| 2 gives the longitudinal beam density profile. Further- 
more, the circular topology of the ring should requires periodic solutions for 'P, with respect to 

x = *(-7rRo,s) and d x <H(nR 0 ,s) = ^(-ttRo, «))• In these conditions, from (7) the 

norm squared (A/ 2 ) of bwf, defined as 


r+xRo 

^ = j_ rRo |^(a:,s)| 2 dx , (g) 

is conserved ( U is assumed a real function), and it has been fixed for simplicity equal to Rq. This 
result is compatible with the physical requirement that A (x,s) dx = N. However, in the 

following we restrict our analysis to consider bunched beam whose effective length is much smaller 
than Rq. Under this assumption, the above conditions of periodicity of do not have a relevant 
role, because in this limit, for the bunch, the ring looks like an infinite linear accelerator. Thus 
we can define the effective bunch length a L , and the expectation value of the momentum V as 


°l{s) 





r+xRo 

dV 

2 

J—itRq 

dx 

dx 


I 1 / 2 


and in complete analogy to quantum mechanics the uncertainty principle 


°l° p l > e L /2 

holds. Furthermore, note that (7) 1/rj plays the role of an effective mass. 


(9) 


( 10 ) 


2«1 The interaction potential and synchrotron oscillations 

As it has been shown in Ref. [4], the potential U can be split in two parts (RF + Self-interaction) 
and (4) becomes: U = U RF + U s . Note that whereas in general U RF is a known function of x 
and s, the explicit expression for U$, depending on the bunch density (collective effects), needs 
particular assumptions about the beam interaction with the surrounding medium to be specified. 
This interaction can be parametrized in terms of the longitudinal coupling impedance [7]. In the 
special case of a linear approximation for the RF-potential, U RF = l^x 2 , where K is the cavity 

strength {^\K\ is the synchrotron wave number), and for a purely reactive longitudinal coupling 
impedance X, the equation for bwf becomes 


d * 
ULT] -d7 


e 2 L r) 2 d 2< S> 
2 dx 2 


+ ^x 2 * 


qi 



(11) 


where I is the beam current, and n is the so-called harmonic number [8]. In the simplest case of 
(X — 0) and for Rq » a L , the Eq. (11) can be exactly solved, and the normalized solutions for 
bwf are the well known Hermite-Gauss modes as it occurs in complete analogy for electromagnetic 
optics in paraxial approximation [9] 


'^m(*£> ^) — 
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In (12) the functions a L (s), p L (s) and <p L (s) are solutions of the following system of differential 
equations 

#OL 

ds 2 


o 

II 

1 

b 

* 

+ 

(13) 

1 1 do L 

(14) 

pL &L dS 


d<t>i _ (Ln 
ds 4 o\ 

(15) 


and H m (x ) are the Hermite-polynomials with m a non-negative integer. Note that |tf m | 2 for m = 0 
(fundamental mode) gives a Gaussian particle distribution. Remarkably, it can be easily proved 
that (13) is completely equivalent to 


tfirr 2 

— + AKo\ — 4 £ with £ 
ds 2 



eirj 2 1 2 

+ t* + 2 K ° l = C(m9t - 


(16) 


In this form it is easily to recognize that the equation for a L (s) (16), i.e. (13), describes the 
synchrotron oscillations. The equilibrium condition $(J l/ ds 0 gives 



.po _ M* PO 
l ~ ~ a L ’ 


( 17 ) 


where a° L and of are the equilibrium value of a L and o p L , respectively, and v. is the synchrotron 
number given by the ratio between the synchrotron frequency fl t = CyJ\K\ and the revolution 
frequency o> 0 [10]. Equation (17) recovers the well known relationship between the bunch length 
a° L and the energy spread af [10]. Since for the present case the bwf is Gaussian, in obtaining 
(17) we have introduced the minimum value o\oi = cl/ 2 of the product aio v L consistently with 
disequality reported in section 2. 


3 Coherent stability criterion and soliton solution 


In this Section, we develop, within the framework of the thermal wave model, an analysis of some 
collective effects occurring when the bunch interacts with the surrounding medium. To this end, 
we consider the special case of RF cavity off and take into account both the space charge effect 
and a purely inductive coupling impedance. Consequently, (11) becomes 


ie L r) 


3* 

ds 


e\r? 3 2 $ 
2 dx 2 


q j 

r£ 0 



(18) 


Note that (18) is formally identical to the cubic NLS equation which describes the propagation of 
an e.m. pulse through a nonlinear medium in paraxial approximation [11], [12], [13]. In this analogy, 
the factor e L ri plays the role of the diffraction parameter (the inverse of the wave number), s 
corresponds to the time, and -[r?g/(X/n)/(27r£ 0 )]|'lf corresponds to a nonlinear refractive index. 
Thus, an analysis of the bunch coherent instability (stability) can be made in complete analogy 
to the electromagnetic case [11], [12], To this aim, applying the well known method developed in 
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nonlinear e.m. optics to search for the sufficient conditions of modulational instability, we show 
that coherent instability for particle bunches is fully equivalent to modulational instability for e.m. 
bunches. Moreover, a soliton-envelope solution, very interesting for accelerator physics, is found. 

As an example, we analyze the instability of a plane-wave (\I> Q (x,s ) = p 0 exp[i(k 0 x - n 0 s)], 
where po is a positive constant) solution of Eq. (18), when a small perturbation around it is 
introduced. Let 

$(r, s) = [p 0 + Pi(x, s)] exp[i(fc 0 z - flo«) + Wi{x, s)] (19) 

be the perturbed solution, being pi and 9\ real functions, and 




qj 

2tt€ L E 0 



( 20 ) 


In order to obtain the dispersion relation of the system we can assume 


pi(x,s) = p°! cos (kx -tts + 4>o) 0i(x,s) = 0?sin(fcx-fts + 0o) , (21) 

where. Pi , 0®, k, ft and <p 0 are real constants. By imposing that (19) is a solution of the linearized 
(18) we obtain the following dispersion relation 



2A?o ± 


( k 2 W 

\ ml 2vE o 



( 22 ) 


Reminding that unstable modes occur for S(ft) ^ 0, namely 


k 2 < 


_ 2 qi_(X\ 2 

nr]€ 2 L E 0 \ n / 0 


(23) 


we get stability for rjX < 0 and instability for r)X > 0. This result recovers the well-known 
condition for coherent stability (instability) for monochromatic charged particle beams [14], in 
addition we remark that the above condition is fully similar to the Lighthill criterion, valid for 
modulational instability of an e.m. plane-wave travelling in a nonlinear medium [11], [12], [13]. 

For a bunched beam (a L « Rq), a solitary solution of (18) is found by looking for a solution 
of a relativistic envelope form: 


tf(x,s) = G(x - /3 0 s)e ix ° x - iwoS , (24) 

w Rh Xo, wo, and /? 0 real numbers. Thus, according to the general theory of NLS equations [11], 
the following soliton-like solution for the beam density (A = AoG 2 ), which satisfies (8), is possible 
under the condition r}X > 0: 


where 


A(x, s ) 


ivy# i 
4elT 0 E 0 tj 
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2e 2 L ToEo T) 
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4 Conclusions 


In this paper, we have reviewed some results and applications of the Thermal Wave Model, showing 
in particular, how it is possible to give a novel approach to the study of the nonlinear longitudinal 
dynamics of a relativistic particle bunch in circular accelerating machines. Neglecting radiation 
damping and quantum excitation, we have shown that the nonlinear interaction between the bunch 
and the surroundings (potential well and wake fields) is governed by an appropriate NLS equation 
(equation (11)), fully similar to the one that holds for the propagation of an e.m. bunch in a 
nonlinear medium in paraxial approximation [11], [12], [13]. Much remains to be done — like, for 
instance, the extension to 2- or even to 3-D, or the development of an iterable formulation — 
to make this model really interesting for the study of the typical, still unsolved, beam-dynamics 
problems. Nevertheless, its very innovative feature of allowing the treatment of the whole beam 
at the same time, makes it look extremely promising for a new, and more complete approach to 
the subject. 
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